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We derive the time evolution of the density contrast to all orders of perturbation theory, by solving 
the Einstein equation for scale-invariant fluctuations. These fluctuations are represented by an 
infinite series in inverse powers of the radial parameter. In addition to the standard growing modes, 
we find infinitely many more new growing modes for open and closed universes. 
PACS numbers: 04.25.Nx, 98.65.-r, 98.80.-k 



The large scale structure of the universe is believed 
to have grown, due to gravitational instability, from 
small primordial density fluctuations. These fluctua- 
tions arc fully characterized by the density contrast, 
5(t, x) = 5p(t, x)/ p(t). At any given time, the connected 
correlation functions of the density contrast, determine 
the power spectrum and the space distribution of the 
matter in the universe [1—5]. The time dependence of 
the density contrast, on the other hand, determines if 
the inhomogeneities would indeed grow, or would just 
oscillate or decay, in the first place. 

In this article, we are primarily concerned with the 
time evolution of the density contrast. In almost all 
standard methods, this is determined by the continuity, 
Euler and Poisson equations in the Newtonian regime 
[1,6-10] and by the Einstein equation in the relativistic 
era [1,6—10]. In the nonlinear regime, the perturbation 
is carried to second order and new growing modes con- 
taining nonlocal terms are obtained for dust universes 
[1,5,11] . The local nonlinear terms have also been studied 
to all orders for the flat universe [4,12]. A single complete 
treatment covering all possible cases and also extending 
to an arbitrary perturbative order does not exist. 

Our starting point is not the usual equations of fluid 
hydrodynamics but the Einstein equation which natu- 
rally covers the present era too. We expand the Einstein 
equation around the background Friedmann universe by 
scale- invariant fluctuations. These fluctuations are ex- 
pressed as an infinite series in inverse powers of the ra- 
dial parameter. This method, which assumes spherical 
symmetry, allows us to obtain the full expression for the 
growth rate of density inhomogeneities in Newtonian, rel- 
ativistic, zero and non-zero curvatures, linear and local- 
nonlinear regimes in a single treatment. Thus in this 
way, we make major simplifications on current pertur- 
bative methods. In a flat universe, the leading growing 
modes are given by an infinite series in the matter and 
radiation -dominated eras. For the open and closed uni- 
verses we find, in addition to the modes presented in the 
literature [8,11], infinitely many more new growing modes 
which are also given by infinite sums. In the nonlinear 
regime, we also find new growing modes for closed and 
open universes in the extreme limits of small and large 
times. 



Contrary to the standard analysis, where the density 
contrast is often taken to be a random Gaussian field, 
we take it to be a scale-invariant quantity. That the 
standard growing modes are contained in our results, is 
a further solid evidence that the time evolution of the 
density contrast is unaffected by its statistical properties. 

Since, the inverse radial parameter now plays the role 
of the perturbation parameter, the order by which the 
structures at different scales are formed is easily deter- 
mined. Our perturbation scheme implies a bottom-up 
or a hierarchical clustering scenario for the formation of 
structures and thus a universe dominated by cold dark 
matter [8]. 

To account for the inhomogeneity of space-time, we use 
the spherically symmetric metric, 
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where the scale factor R p (t, r) is a function of both time 
and coordinate. This metric does not assume homogene- 
ity and is contained in the Tolman metric for a pressure- 
less universe [13], although not restricted to this era. 
We assume that the inhomogeneities are given by scale- 
invariant fluctuations and expand the scale factor as 



R p (t,r)=R(t) + J2 



SR n (t) 



(2) 



around the background Friedmann-Robertson- Walker 
scale factor, R(t). This is a genuine perturbation of 
the metric and not just a gauge mode. A solution £ M 
of 5g^ v = Cfi,u + Cv,n where 5g^ v is the metric perturba- 
tion, such that Sgor = and Sg rr ~ t a r~ n , for a general 
a and n, cannot be constructed [14]. 

The time-time and radial-radial components of the 
Einstein equation are expanded as in (2) and, on using 
the equation of state P — up, are reduced to the follow- 
ing second-order recursive differential equation [15]: 
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5R n 5R m r-( n+m ) = 0, (3) 

where summations over the indices n,m,l and s are im- 
plied and 5R n represents R + 5Rq when n = 0. We have 
solved the above equation by a Maple program to very 
large orders for a flat universe, i.e. k = [16]. The 
solutions are substituted back in the time-time compo- 
nent of the Einstein equation and a series expression for 
the density and subsequently the leading-order growing 
modes for the density contrast are found. The fastest 
growing modes can be written in the closed forms 
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in the matter-dominated era, u = 0, and 
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in the radiation-dominated era. We see that in addition 
to the usual growing modes, t 2 / 3 and t in the matter 
and radiation dominated eras respectively, we have in- 
finitely many more growing modes. The higher growing 
modes are the local nonlinear modes [1,12]. The above 
density contrasts also tell us that smaller structures enter 
nonlinear regime, 6 > 1, faster than the larger structures. 
Thus a hierarchical clustering scenario is expected in this 
scheme. 

For the closed and open universes, the recursive equa- 
tion (3) is highly inhomogeneous and cannot be solved. 
However, assuming that different perturbative orders 
evolve independently of each other, we find, in the 
matter-dominated era, the linearized equation 
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where ip is real for an open and is imaginary ip = —iO, 
for a closed universe. The above equation can be solved 
exactly. The solutions can be used in the Einstein equa- 
tion to obtain the density. The growing modes of the 
density contrast are given by 
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We see that in addition to the usual growing mode [8] 
which is given by n = terms in the above equation, 



infinitely many more growing modes exist in the linear 
regime. The n ^ modes have been overlooked in the 
previous analyses and can have important consequences 
for the power spectrum. 

In the nonlinear regime, the full inhomogeneous equa- 
tion (3) for open and closed universes can also be solved 
analytically at low orders and by a fully algebraic Maple 
program at higher orders [15] in the limits of small and 
large times. At small times, the first growing modes 
S + <~ t 2 / 3 and 5 + <~ t in the matter and radiation- 
dominated eras occur at the lowest order in the pertur- 
bation series, i.e. at the r-independent order [15]. At 
higher orders the solution contains oscillatory and or de- 
caying modes only. Therefore, open and closed universes 
behave flat universe at small times in the linear 

regime, as expected. However, this statement is not com- 
pletely true. In a flat universe the first growing modes 
appear at the third, r~ 3 , rather than the first, order in 
the perturbation expansion, as expressed by equations 
(4) and (5). In other words, at small times the inhomo- 
geneitics will enter the nonlinear phase faster in the open 
and closed universes than in the flat universe. Therefore, 
unlike common expectations, structures can grow faster 
in an open universe than in a flat universe. 

At large times, a new growing mode arises in an open 
universe. Taking the asymptotic limit of the hypergeo- 
metric solutions to the Einstein equation [15] we obtain 
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This result indicates that, at very large times, inhomo- 
geneities will only grow in an open universe if the radia- 
tion pressure is non-vanishing. However, since the bary- 
onic matter cannot grow in the radiation-dominated era 
due to its strong coupling to the radiation, the above 
mode is only relevant for the growth of the perturbation 
in the non-baryonic dark matter component. 

Nonlinear-nonlocal growing modes [1] which would 
arise if radial coordinates, in expression (3), were taken to 
be at different points, shall be discussed in future works. 
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